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DISCREPANCY OF GENERALIZED LS'-SEQUENCES 


MARIA RITA lACO AND VOLKER ZIEGLER 


Abstract. The LS'-sequences are a parametric family of sequences of points 
in the unit interval. They were introduced by Carbone [3], who also proved 
that under an appropriate choice of the parameters L and 5, such sequences are 
low-discrepancy. The aim of the present paper is to provide explicit constants 
in the bounds of the discrepancy of LS-sequences. Further, we generalize the 
construction of Carbone [4] and construct a new class of sequences of points 
in the unit interval, the generalized L5-sequences. 


1. Introduction 


A sequence (xn)neN of points in [0,1) is called uniformly distributed modulo 1 
(u.d. mod 1) if 

1 ^ 

TV ^ 

n—1 

for all intervals [a,b) C [0,1). A further characterization of u.d. is given by the 
following well-known result of Weyl [14]: a sequence (a;„)„gN of points in [0,1) is 
u.d. mod 1 if and only if for every continuous function / on [0,1) the relation 

1 ^ r 

holds. 

An important quantity called discrepancy is introduced when dealing with u.d. 
sequences. It measures the maximal deviation between the empirical distribution 
of a sequence and the uniform distribution. Let lon = {a^i, ■ • ■, a:Ar} be a hnite set 
of real numbers in [0,1]. The quantity 


Dn{ujn) = sup 

0<a<6<l 


N 

ivS 

n—1 


^ > l[a,6)(a;„) - (6- a) 


is called the discrepancy of the given set lon- This definition naturally extends to 
infinite sequences (Xn)neN by associating to it the sequence of positive real numbers 
d^N({xi, X 2 ,..., X]v}). We denote by Dj^{xn) the discrepancy of the initial segment 
{xi,X 2 ,..., xn} of the infinite sequence. 
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Sometimes it is also useful to restrict the family of intervals considered in the 
definition of discrepancy to intervals of the form [0,a) with 0 < a < 1. This leads 
to the following definition of star-discrepancy 

D%{ujn) = sup 

0<a<l 

It is a well known result that a sequence of points in [0,1] is u.d. if and 

only if 

lim D*^{xn) = 0 . 

N—^OO 

Sequences whose discrepancy is of order 0{N~^ logiV) are called low-discrepancy 
sequences. These sequences are of particular interest in the theory of numerical 
integration and are used in the Quasi-Monte Carlo (QMC) integration. For more 
information on discrepancy theory, low-discrepancy sequences and QMC integration 

see [7l[Tl]. 

We now recall the splitting procedure that gives rise to the generalized LS- 
sequences as sequences of points associated to a sequence of partitions of the unit 
interval. 

Definition 1 (Kakutani splitting procedure). If a G (0,1) and tt = '■ 1 < 

i < k} is any partition of [0,1), then air denotes its so-called a-refinement, which 
is obtained by subdividing all intervals of tt having maximal length into two parts, 
proportional to a and 1 — a, respectively. The so-called Kakutani’s sequence of par¬ 
titions (a"a;)„GN obtained as the successive a-refinement of the trivial partition 
^ = {[ 0 , 1 )}. 

The notion of a-refinements can be generalized in a natural way to so-called 
p-refinements. 

Definition 2 (p-refinement). Let p denote a non-trivial finite partition of [0,1). 
Then the p-refinement of a partition tt of [0,1), denoted by pir, is given by subdi¬ 
viding all intervals of maximal length positively homothetically to p. Note that the 
a-refinement is a special case with p = {[0,a), [a, 1)}. 

By a classical result due to Kakutani m, for any a the sequence of partitions 
(Q!"w)„gN is uniformly distributed, which means that for every interval [a, b] C [0,1], 

^ k(n) 

TTN (O = ^ - a, 

n->-oo Kin) ^^ 

^ ^ i—1 

where k{n) denotes the number of intervals in a'^uj = {[t"_i,t"), 1 < * < k{n)}. 
The same result holds for any sequence of p-refinements of w, due to a result of 
Volcic [13] (see also [HIS]). 

The generalized LS'-sequence of partitions represent a special case of a p-refine- 
ment. 

Definition 3 (Generalized LS'-sequence of partitions). Let Li,...,Lk be non¬ 
negative integers, with LiLk 0. We define the generalized LS-sequence of parti¬ 
tions (p2j^ ^^uj)nen BIS the successive p-refinement of the trivial partition to, where 
PLi,...,Lk consists of Li -\- L 2 -\- ■ ■ ■ -\- Lk intervals such that Li has length j3, L 2 has 
length ff' and so on up to Lk having length . 


1 ^ 

^^l[o,a)(a:„)-a 
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Note that necessarily Li/3 + • • • + Lfc/3^ = 1 holds, and consequently for each 
fc-tuple (Li,..., Lk) of parameters satisfying the assumptions made in Definition [3] 
there exists exactly one positive real number /3 satisfying Lij3 + • • • + LkP^ = 1. 

If fc = 2, then we obtain the definition of the classical L5'-sequence of partitions 
introduced by Carbone [4] . Given a sequence of partitions we can assign a sequence 
of points by ordering the left endpoints of the intervals in the partition. The 
corresponding L^-sequences of points, denoted by have been introduced 

by Carbone [4], who proved that whenever L > S there exists a positive constant 
ki such that 

(1) ■ • ■ > ^ ^ 

One purpose of the present paper is to give an estimate for ki. This problem 
still open in [4] can be solved by the numeration approach presented by Aistleitner 
et.al. [2]. However the main focus of the present paper is to generalize the construc¬ 
tion of classical LS'-sequences of points by Carbone [4]. Therefore we introduce a 
numeration system for the integers which is the generalization of the numeration 
system introduced by Aistleitner et.al. [5]. This numeration system will lead us 
to the definition of generalized LS-sequences of points. This definition by a nu¬ 
meration system allows us to prove several auxiliary results on the distribution of 
generalized LS'-sequences which were proved in the case that fc = 2 by Carbone [4] 
and Aistleitner et.al. [2]. All this is the content of the next section. In Section [3] 
we explicitly compute the constant fci in o and obtain explicit upper bounds for 
the discrepancy of the classical LS'-sequences (see Theorem [T]). The computation 
of the discrepancy in the classical case, i.e. the case that fc = 2 will pave the way 
to compute explicit bounds for the discrepancy of generalized LS-sequences, which 
we compute in the final section of this paper. 

2. Generalized LS-sequences 

Let us consider the generalized LS-sequence of partitions. Then the partition 
consists of intervals having lengths ,5",... and this fact makes 

the analysis of the generalized LS-sequences more complicated, compared to the 
analysis of the classical ones, where only two lengths are considered. We denote by 
tn the total number of intervals of correspondingly by ln,i, • • ■, ln,k 

the number of intervals of the n-th partition having length ,8",..., respec¬ 

tively. 

From a general point of view the only canonical restrictions to the fc-tuple 
(Li,...,Lfc) are that the Li are non-negative integers for all f = l,...,fc such 
that LiLfc ^ 0. However, we are interested in low-discrepancy sequences and we 
will see (Remark [2]) that we obtain low discrepancy sequences if and only if all roots 
but one (counted with multiplicity) of the polynomial -I- • • • -I- LiX — 1 have 

absolute value smaller than one. Therefore we will assume from now on that the 
polynomial L^X^ H— • -|- LiX — 1 has no double zeros and that there is exactly one 
root 0 < ,8 < 1 and all other roots have absolute value less than I. Let us remark 
that excluding multiple zeros is mainly to avoid technical issues and similar results 
would be obtained if we only assume that the unique root 0 < /3 < 1 is simple. 

Furthermore let us note that, if we assume that 


(2) 


Li > L 2 > • • • > Lfc > 0, 
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then LkX'^ + • • • + LiX — 1 is irreducible and is the minimal polynomial of the 
reciprocal of a Pisot-number (see i), i.e. a Ic-tuple (Li ,... ,Lk) which satisfies (l2|) 
also meets our assumptions made above. 

Let us write (3i = P and let P 2 , ■ ■ ■, Pk be the other roots of LkX^ + - ■ ■ + LiX — 1. 
When the coefficients of the polynomial fullfill condition ([^l , then the P 2 , ■ ■ ■, Pk 
are the Galois conjugates of P in some order. Note that since 1//3 is a Pisot-number 
by our assumptions we have that \Pi\ > 1 for i = 2,... ,k. 

In a first step we consider the quantities and ln,i for n > 0 and i = 1,... ,k. 
First we observe that in the n-th partition step the longest ln- 1,1 intervals are 
divided into Li H— ■ + Lj. intervals, where Li intervals have length /3", L 2 intervals 
have length and so on. Hence we have that 

tn = ^n-l,l(Ll Lfe) -|- ln- 1,2 ln-l,k 

= tn-1 + + ■ ■ ■ + Lk — 1), 

^n,l — 1,2 d” 

( 3 ) 


1 — ^n—l,k d“ ^k—Pn—l,li 

for all n > 0. Of course the interval [0,1) yields the initial conditions to = Zo,i = 1 
and Zo ,2 = • • • = lo.fc = 0 and therefore we can recursively compute the quantities 
tm ln,i, ■ ■ ■, ln,k for all n > 0. Since from the recursion point of view the sequence 
{tn)n>o is closely related to the sequences {ln,i)n>o, we define ln,o = tn in order to 
state several of our results in a compact way. Moreover we put Inj = 0 if j > fc or 
n < 0. For our purposes we desire an explicit formula for the quantities ln,i- 


Lemma 1. The sequences {ln,i)n>o, with i = 0,... ,k satisfy the recursion 


“t“ ‘ ‘ ' “1“ Tkln—k,ii k. 

In particular, there exist explicit computable constants Xj^i for 0 < i < k and 
I PI j ^ k such that 


( 4 ) 


k 




Proof. First, we observe that by inserting the last line of ([3]) into the second to last 
line we obtain 

ln,k—l — Lkln—2,1 d" I'k—lln—l,! 

Now inserting this identity into the third to last line of ([3|) and going on we end up 
with 

ln,l = Llln-1,1 d- ■ • • d- Lkln-k,l, 

hence we proved the first part of the lemma for i = 1. Since the characteristic 
polynomial of this recursion is 

- Lk 


there exist constants A^.i, with j = 1,..., A: such that 

k 

ln,l=Y.X,,lPr. 

4 = 1 
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Note that the l//3j for j = 1,..., fc are the roots of the characteristic polynomial. 

Since ln,i = ln+i,i—Liln,i we obtain (U), with = {f3~^—Li)Xj^i for z = 2,..., fc 
and j = 1,..., fc. And since = ln,o = in,i + • • • + ln,k we see that 0 also holds 
for z = 0, with Ayo = Ayi + • • • + Ayfe for j = 1,..., k. On the other hand every 
sequence of the form ([3]) fulfills a recursion of the required form and the proof of 
the lemma is complete. □ 

Let us note that the explicit computation of the A’s is easy for a given fc-tuple 
{Li,...,Lk). Indeed one can compute the values of ln,i for all z = 0, 
and n = 0,1,..., k — 1 by using the recursion (jS]). Therefore 0 gives for each 
i = 0,1,... ,k a linear inhomogeneous system with unknowns Solving for the 
A’s by Cramer’s rule we obtain 
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PA 

PI 


/3r' ... 



PA 

Pk^ 


... 

.. 1 1 


pA 

pA 

[0) Aj^i 

/?? .. 


P° 

PA 

Pi 


.. 

PA 

PJ' 

PA 

Pk" 



pA 

pA 

PA ... 

pA 

Let us also state 

another property of the sequences {In 



which we need at several places in the construction of generalized LS-sequences: 

Lemma 2. We have ln,i + • • • + ln,m > + • • • + ln-i,m+i for all m > 1. 

Proof. 11 m > k this is clear since the sequences ln,m sue all strictly monotone 
increasing with n for m = 1,... ,k and ln,m = 0 li m > k. In case that m = 
1,... ,k — 1 we have 

^z-1,1 + . . . + li-i^m+l 

= U-l,l + {h,l — Lili-i^i) + • • • + — Lmli-l.l) 

’Pi li,m.- 

□ 

Our next step is to introduce the numeration system which will be the basis 
for our construction of the generalized LS'-sequence. Let iV > 0 be a fixed integer 
and choose rz such that tn < N < tn+i. We construct finite sequences {em)o<m<n, 
{r]m)o<m<n, (A^m)o<m<n and (Tm)o<m<n recursively in the following way: 

• First we put = N,Tn = e„ = 1 and r]n = [(A^ - . 

Assume that we have computed the quantities Ni,Ti, and zyz for all indices 

n > i > m > 0. 

• We denote by j the unique integer such that Cm+i = ■ ■ ■ = Cm+j = 0 and 
em+j-ei 7 ^ Oj ia particular if Cm+i ^ 0 we put j = 0. 

• Now we compute and Tm'. 

Pm — “t“ ■ * ‘ “t“ lm,j+2- 
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• If Nm < Tm we put em = rjm = 0- Otherwise we put = 1 and 
r]m = [{Nm - Tm)/lm,l\- 

With this definition we obtain a representation for N of the form 

n 

(6) N = ^{ciTi +r}ili^i). 

i=0 

Obviously G {0,1} for all i = 0,1,...,n. Let us note that we also have 
0 < ryi < Li + ■ ■ ■ + Lfc — 2 for alH = 0,1,..., n. Indeed assume to the contrary 

that r]m > Li -\ -+ - 1 and that Sm+i = ■ ■ ■ = em+j = 0 but €m+j+i ^ 0. 

Then we would obtain in case that j > 0 

Tm+l > Nm+l > Tm + (Ll + ■ ■ ■ + Lfe — l)lm,l 

— lm,l + lm,2 + • • • + lm+j+2 + [Ll + • • • + Lfc — l)/m,l 

^ {Lllm,l + I m,2) ‘ + j+ 2 ) 

— ^m+1,1 T ■ ■ ■ “t“ ^m+lj+l — Tm +1 

a contradiction. In case that j = 0 we similarly have 
Tm+l ^ ^m+1,1 “f ^m+1,2 ^ ^m+1,1 

> Nm > Tm + {Ll + ■ ■ ■ + Lfe — l)lm,l 

= lm,l + ^m,2 + • • • + lm+j+2 + {Ll + • • • + Lfc — l)^rn,l 

> {Lllm,l + lm,2) + • • • + {Lj^ilm,! + lm,j+2) 

— lm+l+ “t“ ■ ■ * “t“ lm+lJ+1 — Tm+l 

again a contradiction. 

The following lemma gives a bijection between the integers and digit-expansions 
of the form given above. Note that the following lemma is a generalization of a 
result due to Aistleitner et.al. for the classical LiS-sequences [21 Lemma 3]. 


Lemma 3. There is a bijection between positive integers and finite sequences of 
the form 

T = ((e„,77„),...,(eo,77o)) 

such that Ci G {0,1}, e„ = 1, 0 < pi < Li + ■ ■ ■ + Lk — 2 for all 0 < i < n, Ci = 0 
implies rji = 0 and for all 1 < m < k — 1 we have that rji > Li + ■ ■ ■ + Lm — 1 
implies Ci+m = 0. 

This bijection is given by 

n 

2=0 


and its inverse 

<^{N) = ((e„, ? 7 „),..., (eo, ryo)), 

where the Ti, Ci and pi are computed by the algorithm described above. 


Proof. In order to prove bijectivity we have to show that for every integer N and 
every finite sequence V we have d>(<i)(A^)) = N and <I>(iI'(2?)) = V. The first equa¬ 
tion is evident from the presented algorithm, i.e. $ is injective. Thus we are left 
to prove that <i)(iI'(T>)) = T>, i.e. $ is surjective. The proof is technical and we 
proceed in several steps: 
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Step I: We show that 

$(iV) = ((e„,77„),...,(eo,??o)) 

such that the and rji for all j > 0 satisfy the conditions of the lemma. Thus we 
prove that $(iV) is well defined. 

Note that from the algorithm it is evident that S {0,1}, e^, = 1 and that 

ei — 0 implies rji = 0. Moreover, in the discussion after ([6]) we have shown that 

0 < rji < Li + ■ ■ ■ + Lk — 2 for all 0 < z < n. Therefore we have to show that 

r]i> Li-\ -+ Lm - 1 implies ei+m = 0. 

We proceed by induction on m. We start with the induction basis m = 1, i.e. 
we assume that rji > Li — 1. Let us assume for the moment that e^+i ^ 0. Then 
we get that 

fVi > (.{Ti + rjili^i > li^i + li^2 + [Li — ~ ^j+1,1 > 

a contradiction and therefore we conclude that e^+i = 0, i.e. we have proved the 

induction basis. 

Now, let us assume that rji > Li + • • • + Lm-i — 1 implies ti+m-i = 0 for all 
1 < m < M — 1 and assume that rji > Li + • • • + Lm — 1. We aim to show that 
Cj+M ^ 0 yields a contradiction. By the induction basis, see the paragraph above, 
we may assume that M >2. Moreover by induction we may assume that ei+m = 0 
for all 1 < m < M. Assuming e^+M 7^ 0 implies Ti = li i + • • • + h^M+i and 
Ti+i = k+i^i + • • • + li+i,M- Therefore we deduce 

Ti+i >Ni+i = Ni > Ti + rjili,i 

> h,i + ■ • • + li,M+i + {Li + ■ ■ • + Lm — ^)h,i 

= {Lili^l + li^2) + • • • + {LmU,i + h,M+r) 

= ^i+1,1 + • • • + li+l,M = Ti+i 

a contradiction, i.e. e^+M = 0. 

Step II: It is enough to show that there are exactly — 1 sequences of length < n 
satisfying the conditions of the Lemma. 

Indeed, we have already seen that <!> is injective. Therefore we have to show that 
is surjective. In particular, it is enough to prove that ih induces a surjective map 
between the positive integers < and sequences of length < n, which satisfy the 
restrictions of the Lemma. Therefore we have to prove the following claim: 

Claim 1. There are exactly — 1 sequences of length < n satisfying the conditions 
of the Lemma. Moreover, there are exactly ln-i,i such sequences of length n of the 
formV = ((1,0),...). 

We will prove that claim be induction. 

Step III: The induction basis is evident, since there are exactly ti — I = Li + 
• • • + Lfe — 1 sequences of length 1. Moreover there is only one sequence of length 
1 starting with the pair (1, 0). 

Step IV: Suppose the claim is true for all integers M < n + I. We show that it is 
also true for M = n + 1. 
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First, let us show that there are exactly ln,i sequences of length n + 1 starting 
with (1,0). Let 1 < m < k, then there are exactly (Li + ••• + Lm — l)^ra-m,i 
sequences of the form 

m—1 times 

((1,0), (0,0),.?., (0,0), (1,7?),...) 

by induction hypothesis and the assumptions of the lemma Similarly there are 
exactly tn-k sequences of the form 

at least k times 

((1,0), (0,0),.^.,(0,0)', (1,7?),...). 

Therefore the number of sequences of length tt, + 1 starting with (1,0) is 

(Ll — l)ln-l,l + {Li + L 2 — l)ln-2,l + . •. 

+ (Ll + ■ ■ • + Lk ~ l)^rt—fc,l + tn—k = 

(Ll — + (Li + L 2 — l)ln-2,l + . . . 

+ {Ll + ■ ■ • + Lfc-1 ~ l)^ra-fc+l,l + tn-k+1 = 

(Ll — 1)Z„_1_1 + (Ll + L 2 — l)ln-2,l + . • . 

+ (Li + • • • + Lfc_i — l)ln-k+l.l + in-k+1,1 + • ■ ■ + ln-k+l,k- 

Now note that 

(Li+ • • • + Lm-l — l)^n,l + in,l + ’ ’ • + ln,m 

= + ^n,2) + ’ ’ ’ + + ln,m) 

= ^n+1,1 + • • ■ + ^ra+l,m-l 

applying this identiy to the equation above we obtain that there are 

(Ll — + (Ll + L 2 — l)/ri-2a + . . . 

+ (Li + • • • + Lfe_i — l)^„_fc_|_i^i + ^„_fc_|_i_i + • • • + In-k+l^k = 

(Ll — 1)^„_1^1 + (Ll + L 2 — l)/n-2,l + . . • 

+ (Ll + • • • + Lk-2 — l)^n-fc+l,l + ln-k+2,1 + ’ • ’ + ln-k+2,k-l = 


(Ll — l)^„_i^i + Z„_l,l + ln-1,2 — Llln-IA + ln-1,2 — in,l 

sequences of length n + 1 starting with (0,1). Therefore there are (Li H— ■ + Lk — 

l)ln,i sequences of lenght n + \ starting with (1, 7 ?) and 0<7?<LiH - + Lk — 2 

and by induction there are — 1 sequences of length < n +1. Therefore all togehter 
there are 

(Ll + ■ ■ ■ + Lk — l)^n,l + tn ~ 1 = tn+l ~ 1 
sequences of length < tt + 1, satisfying the conditions of the Lemma. □ 

The numeration system ([6|) allows us to generate the generalized LS-sequences 
of points in a direct way, as clarified in the following Definition 
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Definition 4. Let N be an integer with representation given in ([6|). Then 


min{Li, e, + f],} 

i^O 

+ /3*+^(max{ei + iji - Li,0} - max{ei + pi - Li - ^2,0}) 

+ /3*+^(max{ei + pi - Li - ^2,0} - max{ei + rji - Li - L 2 - ^3,0}) + 

+ /3*+''“^(max{ei + pi - Li -i/c- 2 , 0} 

- max{ei + pi - Li - Lk-i,0}) 

+/3*+''(max{ei+ ryj - Li-Lfe_i,0})) . 

The easiest example of a generalized LS'-sequence is obtained by considering 
only three possible lengths for the intervals determining the partition. We call such 
sequences LMS-seqnences. 

Example 1. Consider the LMS-sequence of partitions (pf gUj)nen defined by 
the equation /3^ + /3^ + 2/3 — 1 = 0. In the following we describe the sequence of 
partitions and the associated sequence of points obtained by the procedure introduced 
above. 


0 


2/3 2/3 + 1 


++ 


++ 


H 


2/3 


i I /3 
2/3= +/3= 


/3 + /3= /3 + 2;S= j2^ 2/3 + /3= 3 

/3 + 2/3= + /3^ 


The sequence of partitions is obtained by splitting at each step the longest interval 
into 2 long intervals followed by an interval of medium length and by a short interval. 
So at the first partition we split the unit interval into two intervals of length fi, 
followed by one of length /3^ and one of length (3^. At the second step we split only 
the first two intervals proportionally into two intervals of length ftS, one of lenght 
and one of length fi'^, respectively. The procedure goes on in this way as for the 
classical LS-sequences. 

The associated sequence of points is the sequence of left endpoints of the in¬ 
tervals determining the sequence of partitions, where the order is determined by 
Definition\^ 

Let us list the digit expansion of the first 10 positive integers and the correspond¬ 
ing points of the sequence. 
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$(!) = 

((1,0)) 

-)• 

/3 

$(2) = 

((1,1)) 

-)■ 

2/3 

$(3) = 

((1,2)) 

-)■ 

2/3 + ^2 

$(4) = 

((1,0), (0,0)) 

-)• 

/32 

$(5) = 

((1,0), (1,0)) 


/3 + /9" 

$(6) = 

((1,1), (0,0)) 


2^2 

$( 7 ) = 

((i,i),(i,o)) 


/3 + 2/32 

$(8) = 

((1,2), (0,0)) 


2/32 + ^3 

$(9) = 

((1,2), (1,0)) 


/3 + 2/32 + /3' 

4>(10) = 

((1,0), (0,0), (0,0)) 

-)■ 

/33 


Remark 1. Let us point out the connection of the generalized LS'-sequences to the 
classical LS'-sequences and van der Corput sequences: 

The case k = 1: In the case A: = 1 we have /3 = 1/Li and the numeration 
introduced above is the usual Li-adic numeration. Indeed, let us note that 
In = tn = Tn = L" and in particular we obtain 

n n 

N = 'y^{eiTi + diLi, 

i—0 i—0 

where di = ei + rji G {0,1, Li — 1}. A close look on Definition U reveals 
that in this case the generalized LS'-sequence coincides with the van der 
Corput sequence. 

The case k = 2: In this case the generalized LS'-sequence coincides with the 
classical LS'-sequence. Let us note that in the case k = 2 we have that 

^n,2 

n n 

N = 

i—0 i—0 

which correspondes to the numeration system introduced by Aistleitner 

et.al. [2]. 

In order to give an estimate for the discrepancy, it is necessary to introduce the 
notion of elementary intervals. An interval is called elementary if it is an element 
of fo'' some n. Equivalently we can define elementary intervals as all 

intervals of the form some m, where 

$(a; + = ((1, , (eo, %)) 

with rjm-i < Li. In particular, there exists an integer y < tm such that + 

/S'" = Obviously we may choose y = x + lm-i,i < tm- 

The next step consists in finding a method to decide whether a point is 

contained in some given elementary interval or not. 

Lemma 4. Let + /^™) elementary interval. Then 

m— 1 

^ ^ ^ {^i^i “ 1 “ 

i^O 

is the truncated representation of N. 
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In addition let = ^{l : I < iV, G and assume that 

n 

TV = a; + ^ {eiTi + 

i—m 

Then 

n—m 

A^r\N) = Y, {e,+mT, + Tl,+mk,l) + 1 . 

Proof. As soon as we have shown that all integers N of the form 

n 

N = x + Y + r|^h,l) 

i—m 

have the same first m digits as $(x) we can proceed as in the proof of [S] Lemma 4]. 
But, that all those N have the same first m digits is evident from the constraints 
on to be an elementary interval, i.e. x is such that 

$(a; + = ((1, Vm-i), ■■■, (eo, Vo)) 

with Vm-i <471 — 1, and Lemma [31 □ 


Next we are interested in an accurate formula for —where G . 

Lemma 5. Let us assume that N has a representation of the form and assume 


( 7 ) 

where 


(m) 


Then we have 

A)r\N) 

N 


D 

= + — 

^ + AT ’ 


i? < 1 + |Ai,o| + ^(2Aj + |Aj,o|)- 
i=2 


and 


A,- = max ■ 
e=2,...,k 


Ei=i + (4/1 -I- 1 Lk- 2)|Ajp| 

1-1/9.1-' 


Remind that we assume that the polynomial LkX^ + ■ ■ ■ + LiX— 1 has no double 
roots and has only one positive root /3 < 1. 


Proof. Using our assumptions and Lemma |3| we can calculate the exact values of 
aY’{N) and N. In fact, we have 

n n 

N = x+ ''Y. {^jTj + Vih,i) and A^Yi^) — 'Y. {ijTi-m + Vih-m,i) + 1) 

z=m i—m 

where = k-m,i + • — I- k-m,j such that j is the integer such that e^+i = • • • = 
Ci+j = 0 and Ci+j+i 0. In particular, if we write 

k 

T, = ^/3-*A,„ 
i=i 
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then we have 


f. ■ ■ 

i=i 

where Aj^i = A^-.i + • • • + for some integer 2 < £i < k depending on i. This 
yields 

^ + Vik,l) + X 

_ + ’7*^1,l)/?”*''''" + J2’j=2 ^3 + 1 

E^=n^(e^Al.. + 77 *Ai.i)/3-*+- + E-=2 + *’ 


where 


4 = ^(e.A,-, + ,7,A,,i)/3-*+™<A, 

i—m 

for j = 2,..., k. Further, note that 

n k 

N = r^Y. ^ + a; 

i—m i=2 

Therefore we obtain 

A^r\N) , E-=2^i+ 1 , E”=™(e.Ai.. + ^.Ai.i)r^+--fV/3^ 


N 


=/3- + 


TV 


TV 


=/3™ + 


E,t2^i+1 , Er=™(e.Ai., + i7.Ai.i)r*+" 


TV 


TV 


Z]r=m(^iAl,i + 77iAl,l)/3 + Y!j=2 ^3 (;^) + 2^/^” 


TV 


=/3™ + 


1 + E,t2^i(i-(|)™)-^^" 

N 


Finally we want to estimate a;/?"*. Since we assume that ^{x + lm-i,i) = 
((1, rjm-i), ■ • ■, (eo, ?7o)), with 77^-1 < Li — 1 we obviously have that 


X < im —1 H“ ^ 771 — 1,1 (Z/1 1) < tm, 

/3 


and therefore 

k n fn k 

x/3™ < [Ai^ol + |Aj,o| — < |Ayo|- 

T =2 3=1 

If we put all our results together and note that 11 — (/3//3j)™ | < 2 for each positive 
integer m we obtain the statement of the Lemma. □ 


Remark 2. We note that the only place, where we used the fact that the polynomial 

LkX^ H-h LiX — 1 has no double zeros and all roots but one have absolute value 

smaller than 1 is in the proof of Lemma[5l Let us note that dropping the assumption 
that there exists no double zeros would not change the result. But, it would result 
in a slightly worse estimate for R and more technical difficulties in the course of 
the proof. For the sake of simplicity we stick with the case of simple zeros. 
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Further let us note that assuming that all roots but one have absolute value < 1 
would result in an estimate of the form 


\R\ < (logN)^CL,,...,L, 


where A is the number of roots with absolute value = 1 counted with multiplicities 
and is some constant depending on the A:-tuple Plugging 

this bound into the proof of Theorem [2] we would obtain 




N 




J < Cli,...,L2 


(log N) 


A+l 


N 


In case that LkX^ + • • • + LiX — 1 has more than one root with absolute value 
> 1 the bound of the discrepancy would be only of the order N~^ for some B <1. 
Let us note that similar observations have been made by Carbone [4] in the case 
that k = 2. 


In the next section we compute the discrepancy of classical LS-sequences. To 
do so we need the following Lemma proved in [2 Lemma 6] which is more precise 
than Lemma m 


Lemma 6. Assume that N has a representation of the form ([H]), and assume that 
■ we have 


■ T'hen we have 


( 8 ) 

where 

with 


N 


N 


^(em+i?zAi)(-5/3r 


-L-2S + V 1 JTIS , -L + VL^ + 45 

B ^ ^ - Ai = 


2VL2 + 45' 
Moreover R can he estimated by 


2VL^ + 45' 


1 _ ('Cfl'in-m+l 

\R\ < max{|ri|, \r, + {L + S - 2)Ai|}--. 

if SP ^ 1 and 

|i?| < max{|ri|, |ti + (L + S' — 2)Ai |} max{n — to + 1,0} 

ifSP=l. 


3. Discrepancy bounds for classical LS-sequences 

This section is devoted to the explicit computation of the discrepancy of classical 
LS-sequences with the aim proving 

Theorem 1. Let 5 )neN a classical-LS-sequence of points with L > S. Then 


(9) 


Rn{^L,S^ ^L,Sy • ■ ■ ! ^L,s) ^ 


log 

N\ log^l 


(2L + S - 2) 


R 

i-sp 



B 

TV’ 
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with R = max{|Ti|, |ti + (L + S' — 2)Ai|} and 


B = {2L + S-2) 


R 

l-S/3 



+ 2 


Proof. Let Ci Ci s; ■ • ■; s ^ points of the sequence (^£ s)nGN and 

take n = max{fc : tfe < iV < tk+i}- Consider an arbitrary subinterval [x, y) in [0,1[. 
We want to estimate the number of points among Ci sj s; ■ • ■; s "''^hich lie in 
the interval. To do so we approximate the interval [x, y) from above by elementary 
intervals of length at least and apply to each interval. The points x and y 
belong to some interval determining the n-th partition respectively. In particular let 
us assume that x lies in an elementary interval [x",x') and y lies in an elementary 
interval of length /3" respectively. 

We approximate the interval [x, y) from above, i.e. we try to cover the interval 
[x'^ y') C [x, y) by as few as possible, disjoint elementary intervals of length at least 
jjri+i_ Starting from the point x' we move to the right in order to reach the point 
y' with steps of variable length. By the definition of our LS'-sequence there is an 
integer Q < ^ < L — 1 such that x'l = x' + £/3" + 5'/3”+^ is the left endpoint of an 
elementary interval of length at least In case that y' < x\ we have 

found our covering. On the other hand there exists an integer 0 < ^ < L — 1 such 
that y\ =y' — ifdP is the right endpoint of an elementary interval of length at least 
Hence we are left by the problem to cover the interval {x'i,y'i) by as few as 
possible, disjoint elementary intervals of length at least /S". It is now easy to see 
by induction on n that [x', y') is covered by at most 2L + S — 2 disjoint elementary 
intervals of length /3^,... ,/3" respectively and at most S intervals of length 
and L invervals of length /3. 

Let Ax,y{N) = jl{Z : I < ^ [x,y)}- Since by construction there are no 

[x', x) U [y, y') with I < N we obtain by Lemma El 


( 10 ) 


(y-x)- 


■^x,y{N) 


N 


<|x - x'l + \y- y'\ + 


{y - x') - 


Ax',y'{.N) 


N 


n+1 


<2/3” + (2L + S' - 2) ^ 


m—0 


i?(l - (-S/?)™) +1 


N 


where we used the inequality 

1 > 1 - x/3™ > 1 - = 1 - Ao - Ai(-S/32)™ > 1 - Ao - Ai = 0, 

where Im = lm,i = Aq/?™ + Ai(—S/S)’" and Aq = order to establish 

Theorem [T] we are left to estimate the sum in (flUl) . Since |i?| < 
find 


n+1 

^|i?(l-(-S/3)+| 

m=0 


n+1 


m—O 


(1 - (-S/3)’")(l - (S/3)”-’"+i) 
1 - (5/3) 


n+1 

<^E 


m=0 


(1 + (S/S)™))! - (S/3)"-’"+i) 
1 - (5/3) 


^ (n + 2)R 

- I-Sj3 
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Since (3 " < we get n < and we obtain Theorem [TJ □ 


As an example we can compute the discrepancy of a particular LS'-sequence 
obtained by taking L = S = 1. This sequence, also called Kakutani-Fibonacci 
sequence, has been also analyzed in detail in UM in the frame of ergodic the¬ 
ory where it has been shown that it can be obtained as the orbit of an ergodic 
transformation. 

Take L = S = 1, then by Theorem [T] we have that 




2.366 


^ogN 

N 


3.139 

N 


In the case that L = 10 and A = 1 we obtain 


Ci0,1j • ■ ■ ! ClO,l) — ^-66 


AogN 

'W' 


22.02 

N 


In particular we obtain 


Corollary 1. Let S be fixed and assume that L is large. Then we obtain 

■ I 2L 


lim 

N —¥00 


logN 


\ogL 


as L 


Proof. By the formula given for R in Lemma [6] and the fact that /3 
we obtain that R ^ SjS ^ S/L hence 


_ -L+VL^+45 
“ 2S 


2L + S ( R \ 2L 

|log/3| yl - 5^ j \ogL 

as L —>• 00 and S is fixed. □ 


Remark 3. One can rather easily improve our bound for the discrepancy by a factor 
1/2. This can be done by balancing our choice of the intervals that cover [x,y). 
Let us assume that x' is the right endpoint of an interval of length /3" and let us 
assume that x' lies nearer to a right endpoint x'l than to the left end point and we 
go from x' to the right instead to the left and proceed in this manner, we need fewer 
intervals (roughly 1/2-times fewer) to cover the interval \x',y'). We did not work 
out the details for this improvement, since the paper is already rather technical and 
this approach would further increase the technical difficulties. 


Remark 4. Fame [8] (see also m page 25]) obtains the upper bound 


lim fVZ5^(^i,...,e^) 

Af —>00 


L 

41ogL 


if L is large, where denotes the van der Corput sequence for base L. This shows 
that our approach gives up to a factor 8 (respectively 4 considering the remark 
above) a similar main term for the discrepancy as one obtains for the van der 
Corput sequence. 


At the end of this section we want to compare our approach to obtain explicit 
bounds for the discrepancy of LS'-sequences to the approach due to Carbone [4]. 
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Therfore we give explicit bounds for the star-discrepancy of the Kakutani-Fibonacci 
sequence of partitions Dn(ppi) and of points where 


Dnipll) 


sup 

0<b<l 


Mo,b)i.ti) - b 

i=l 


Following [4], in order to find upper and lower bounds for Dn{pi i), we need to 
estimate 

. T 

i—l 

Thus we consider [0, h) as a union of intervals defining the p-th partitions p\ for 
p < n—\. Therfore let us count how many consecutive intervals • Imi of Pi,i 

are contained in [0,5). Now, we count how many consecutive intervals /J, 
of p\ i are contained in [0,5) \ Ufc\ so on. Of course it may happen that 

TOfc = 0 for some k. Going on with this procedure we get 


u 

2=1 


/r'c[o,5)\u 



Thus 

n-l /rrip \ 

[o.^)=u ■ 

p=i \i=i j 

In particular, we need to compute rup. Since L = S = 1, we have that if If = 
LP = JP, then /f U If = If~^ for some i, and if /f = = Lp, then If = Sp~^. In 

particular, this shows that rUp = 1. Thus there are no consecutive intervals of the 
same partition. After some tedious computations, as in [4l Equation 10], we finally 
get 


■ T 

2=1 

= - ^ + + /5) - (-/3)”) ■ 

By considering separately the case that n is even and the case that n is odd, one 
gets the following bounds: 


and 


Thus 


0.0652 1 ^ , 0.4068 

—- <—l^I[Q,b){U)-b<— - if n even 

2=1 


■ 


-0.2764 1 ^ , 0.2764 

<-T.Mo,bm-b<^- 


■ 


2=1 


if n odd . 


Dniph) ^ 


C2 


with C 2 = 0.4068. 
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Now, let (Pl,s)»ign be the sequence of long intervals In of p2,s- Then it is possible 
to bound its discrepancy (see [H Proposition 3.4]) and obtain for ^ < L + 1 and 
every u G N 


( 11 ) r ^ < r ’ 

where Ci and 62 are constants independent of n. Furthermore, it can be shown that 
if S' < L + 1, then 

, ^Is) < 2 "' ^ %~logP\ logIV + (L + S - 2) , 

where C 2 is the constant in the upper bound of Dn{p^ 5 ) (see [U Equation 29]). 

Finally, plugging in L = S = 1, we can give an estimate for the discrepancy of 
the sequence (^"i)nGN. 


1^a(0 < 2 


C2 log N 
llog^l N 


2 0.4068 logiV 

|log/3| N 


1.6911 


logiV 

N 


Let us point out that this result is in accordance with the following relation 


D%{Xn) < DN^Xn) < 2D%{Xn) , 
which holds true for every sequence of points (a;„)„gN. 


4. The Discrepancy of generalized LS-sequences 


The aim of the present section is to provide bounds for the generalized-LS- 
sequences: 


Theorem 2. Let ^^)neti be a generalized-LS-sequence of points. Then 

(12) log(iV+l)-log|Ai,o/3'=' 


iV|log/3| 


-(2Li + L 2 + ■ ■ ■ + Lk — 2 )i?, 


where 


and 


R — 1 + |Ai,o| + ^( 2 Aj + |Aj,o|) 
1=2 


A, = max < 
t=2....,k y 

provided N is large enough. 


I + (Ti H- + Lk- 2 )|Aj 4 | 

I-I/37I-' 


Proof. The proof of Theorem[^runs along the same lines as the proof of Theorem[T] 

As above we put n = max{m : tm ^ N < tm+i} and we consider an arbitrary 
subinterval [a:, y) C [0,1) and approximate it from above by elementary intervals. 
As in the LS-case we assume that x G [x'',x') and y G [y',y'') lie in elementary 
intervals of length /3" respectively. 

Let x\ be the next left endpoint of an elementary interval of length at least ,0"“^ 
from x'. We need at most Li + • • • + — 1 intervals of variable length to cover 

the interval [x',Xi). Similarly we proceed for the interval [y(,j/'), where is the 
nearest right endpoint of an elementary interval of length at least lying left 

of y'. Obviously we need at most Li — 1 elementary intervals of length /J" to cover 
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[yi,y'). Similar as in the proof of Theorem [T] we conclude that we need at most 
2 Li + L 2 + • —h Tfc — 2 elementary intervals of each length /3^ with £ = 1,..., n + fc 
to cover [x', y'). 

Using similar notations as in the proof of Theorem[T]and by Lemma[5]we obtain 


(13) 

where 


{y-x)- 




N 


<|a;-a;'| + \y-y'\ + 


(j/'-cr')- 


N 


^2/3^ + (2Z/1 + ■ ■ ■ + -Lfc — 2)(n + /c)_R, 


and 


^ — 1 + |Ai,o| + ^(2Aj + |Aj,o|) 
1=2 


A _ 1 I Z)i=i + (^1 H- 1- Ufe - 2)|Aj,i| 

A, - max < l-|/3,|-i 


e=2,...,k 

Since |Ai,o|/3“" — ^ < tn < N, provided N is large enough, we get 

log(Af + l)-log|Ai,o| 


n < 


llog^l 


and we obtain Theorem [2j 


□ 


Let us compute the discrepancy for a concrete example: 


Example 2. Let us consider the case of the LMS-sequence considered above (see 
Example m, with L = 2 and M = S = 1. To ease the notations we write In = 
Inp^ixin = In ,2 o,nd Sn = In, 3 - First we note that the roots or the polynomial 
X3 _|_ 2X — 1 are (3 ~ 0.393, (32,3 — —0.696 ± il.436. Moreover, we can apply 
Cramer’s rule as in to compute the Xj^i’s, for 1 < j < 3 and 0 < i < 3. Then 
by Theorem\Mwe have that 


^Nif.2,1,1^ ^2,i,ii • ■ • I ?^i,i) — 51.4562 


\og{N + 1) 
N 


122.5173 

N 


Finally let us state the following remark. 


Remark 5. Let us note that the proof of Theorem [5] leaves a lot of room for im¬ 
provement. First, as already explained in Remark [3] by more carefully choosing the 
intervals which cover [x', y') we might replace the factor 2Li -|- L 2 + • —I- Lfe — 2 by 
something like Li -I- L 2 -I- • —h Lfe. Moreover, the estimates for Kj and R in Lemma 
[ 5 ] are rather rough and can be certainly improved in concrete cases. 
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